Giant thermopower and figure of merit in single-molecule devices 
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We present a study of the thermopower S and the dimensionless figure of merit ZT in molecules 
sandwiched between gold electrodes. We show that for molecules with side groups, the shape of the 
transmission coefficient can be dramatically modified by Fano resonances near the Fermi energy, 
which can be tuned to produce huge increases in S and ZT. This shows that molecules exhibiting 
Fano resonances have a high efficiency of thermoelectric cooling which is not present for conventional 
un-gated molecules with only delocalized states along their backbone. 
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Recent advances in single-molecule electronics have 
been underpinned by considerable improvements in elec- 
trical contacting techniques'^, which have allowed the 
identification of a range of fundamental properties includ- 
ing switching-^, rectification^, memory^ and sensing 9 . 
Single-molecule devices offer a potential route to sub- 
10 nm electronics. For the purpose of designing large- 
scale integration of such devices, a knowledge of ther- 
mal as well as electrical properties is needed. Further- 
more, a knowledge of thermoelectric properties of single 
molecules 10 may underpin novel thermal devices such as 
molecular-scale Peltier coolers and provide new insight 
into mechanisms for molecular-scale transport. For ex- 
ample from the sign of the thermopower it is possible 
to deduce the conduction mechanism or relative position 
of the Fermi energyi^, with a positive sign indicating p- 
type conduction, which means the Fermi energy is lying 
closer to the HOMO level. The thermopower can also 
show very interesting behaviors, like a linear increase 
with the molecular length—. Other theoretical studies 
have shown the possibility of changing both the sign and 
magnitude of the thermopower by electrically gating a 
molecular wir o 14 i 15 , which moves the Fermi energy across 
a transmission resonance and thus changes its sign. 

However when the Fermi energy lies within the 
HOMO-LUMO gap, thermopower is expected to be low. 
Indeed recently, Reddy et al. 11 measured the room tem- 
perature thermopower of 1,4-biphenyldithiol, using a 
modified scanning tunnel microscope. They found the 
thermopower to be +12.9 /iVK -1 . This low value was at- 
tributed to the Fermi energy sitting within the molecule's 
HOMO-LUMO gap where both the DOS and transmis- 
sion coefficients are relatively flat. 

In this paper we show that much higher values can be 
obtained from molecules exhibiting Fano resonances. We 
investigate the thermoelectric properties of two different 
molecules, the 1,4-biphenyldithiol (BPDT) and CSW- 
470-bipyridine (CSW-479) 16 , as a function of the rota- 
tion angles of some of the molecular groups. Both ex- 
hibit Breit-Wigner-like resonances in the transmission co- 
efficient T(E) near the Fermi level corresponding to the 
HOMO and LUMO levels. However, CSW-479 also pos- 
sesses states associated with the bi-pyridine side group, 
which are weakly coupled to the backbone and produce 




FIG. 1: Composition of the studied molecules a) 1,4-Biphenyl- 
dithiol (BPDT) and b) CSW-479-bipyridine (CSW-479) 



an additional Fano resonance in the electron transmis- 
sion coefflcient a 17 i 18 . The position of the Fano resonance 
can be tuned via a variety of mechanisms, including ring 
rotation of the side group^ and gating by nearby polar 
molecules such as water 19 , whereas Breit-Wigner reso- 
nances are relatively insensitive to such effects. In what 
follows, we demonstrate that changing the conformation 
of the side group in CSW-479 moves the Fano resonance 
close to the Fermi energy and creates huge changes in the 
magnitudes of the thermopower and can even change its 
sign, whereas only small changes occur when the confor- 
mation the BPDT molecule is changed. 

The thermal and electrical properties of a molecular 
device can be understood by studying the thermoelectric 
coefficients G, L, M and K, which relate the charge cur- 
rent / and heat current Q to the electrical bias AV and 
the temperature difference AT across the system. These 
equations, I = GAV + LAT and Q = MAV + KAT, 
are commonly rewritten in terms of measurable thermo- 
electic coefficients: the electrical conductance G, ther- 
mopower 5, Peltier coefficient n and thermal conduc- 
tance ft, defined by: 
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Thermoelectric properties can therefore be utilized in 
either charge-driven cooling devices or heat-driven cur- 
rent generators. In conventional devices the maximum 
efficiency of either heat transfer or current generation 
is proportional to the dimensionless thermoelectric fig- 
ure of merit ZT — - ^ T defined in terms of the mea- 
surable thermoelectric coefficients, where high efficiency 
corresponds to values ZT » 1. Some of the high- 
est measured and predicted values of S and ZT are 
found in materials that have sharp densities of states 
(DOS) 2Q i 21 such as bismuth nanowires^ 2 - or heavily doped 
semiconductors 2 - 3 . It is therefore expected that similar 
high values could be obtained in molecular devices since 
molecules also have sharp DOS close to their molecular 
levels. 

The thermoelectric coefficients 5, G and k can be com- 
puted by extending the Landauer-Biittiker formalism to 
include both charge and heat currents. In the linear tem- 
perature and bias regime the three required thermoelec- 
tric coefficients 24 can be expressed in terms of the mo- 
ments L n of the transmission coefficient via the following 
equations 25 . 




L n = I°°JE - Ef TT(E)^±dE (6) 

If T(E) is a slowly- varying function of E, then at low 
temperatures, theses expressions can be simplified by ex- 
panding T(E) about E = Ep. In the present case, where 
T(E) varies rapidly when a Fano resonance is present, 
we calculate the full integrals rather than taking the low 
temperature limit. 

To calculate thermoelectric properties from first princi- 
ples, we use the density functional theory code SIESTA 26 
which employs norm-conserving pseudopotentials to re- 
move the core electrons 2 ^ and a local atomic orbital basis 
set. In particular we use a double-zeta polarized basis 
set for all atoms and the local density approximation 
for the exchange and correlation functional as parame- 
terized by Ceperley and Alder 28 . The Hamiltonian and 
overlap matrices are calculated on a real space grid de- 
fined by a plane wave cut-off of 200 Ry. Each molecule 
is relaxed into the optimum geometry until the forces on 
the atoms are smaller than 0.02 eV/A and then sections 
of the molecules are rotated to the desired angles. The 
contact atoms are attached to the hollow site of a gold 
(111) surface consisting of 6 layers with 25 atoms per 
layer. The transmission coefficients are found by using 
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FIG. 2: Transmission coefficient for the a) BPDT and b) 
CSW-479 molecules for twist angles of 0° (black), 30° (red), 
and 60° (blue). 



the non-equilibrium Green's function formalism 29 . 

The large values of S and ZT predicted below rely 
only on the ability to move the Fano resonance close to 
the Fermi energy and are independent of the physical 
mechanism used to produce this shift. Furthermore, even 
though DFT may not predict the exact position of the 
Breit-Wigner and Fano resonances, the fact that the side 
groups are also aromatic groups and have energy levels 
close to the energy levels of the backbone, ensures the 
presence of Fano resonances near the Fermi level. There- 
fore an adjustment to the environment of the side group 
will generate the required shift. In what follows we exam- 
ine the effect of rotating either the phenyl rings in BPDT 
or the side group in CSW-479, which could be induced 
by tilting the molecule using an STM tip 3 or via steric 
hindrance 3 ^^^ 2 -^ 3 -. 

The sections of the molecules to rotate are chosen so 
that in BPDT the coupling along the backbone and in 
CSW-479 the coupling to the side group (see Fig. (pQ)) 
are changed. For BPDT the molecule is first relaxed 
and then the torsion angle between the phenyl rings is 
artificially increased from 0° to 90°, where 0° corresponds 
to the case where the rings are planar. For CSW-479 the 
side group is rotated from 0° to 180° about the C-C bond 
that couples it to the backbone, where 0° corresponds to 
when the bi-pyridine side group is approximately in line 
with the molecular backbone. 

Fig. ([2]) shows the transmission coefficients for a) 
BPDT and b) CSW-479 as a function of energy for dif- 
ferent twist angles when the molecules are connected to 
the hollow site on the gold surfaces. As the phenyl rings 
are rotated from 0° to 60° in BPDT, the conjugation 
of the molecular backbone is reduced, which moves the 
LUMO level up in energy and reduces the transmission at 
the Fermi energy, following a known cos 2 (#) la w 31 i 32 i 33 . 
However the overall shape of the transmission remains 
almost unchanged. In contrast the CSW-479 molecule 
possesses a Fano resonance near the Fermi energy and as 
the side group is rotated, the Fano resonance moves to a 
lower energy, leading to a significant change in the slope 
of T(E) near E = E F . 

For each twist angle, Fig. ([3]) shows the thermopower 
for a) CSW-479 and b) BPDT, obtained by using the 
transmission coefficients shown in Fig. (2), to evaluate 
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FIG. 3: Contour plot of the thermopower S for a) CSW-479 
and b)BPDT as a function of temperature and twist angle. 
Note that the color scales used in the two figures differ by 
orders of magnitude. 




FIG. 4: Contour plot of the figure of merit ZT for a) CSW-479 
and b) BPDT as a function of temperature and twist angle. 

equations ([2]) and ([6]). The thermopower in the BPDT 
system increases smoothly with rotation angle at high 
temperatures and remains positive. This is expected 
since around the Fermi energy the transmission coeffi- 
cient decays slowly for all angles and therefore the asym- 
metric term Li in Eq. (2) is always negative. The 
thermopower for the optimized geometry with ~ 35° 
is calculated to be 6.49 /iVK -1 , which is close to the 
value measured in Ref. 11. The discrepancy could be at- 
tributed to differences in the contact geometry. In fact, 
we know from previous calculations that by altering the 
contact configuration or by varying the distance between 
the leads and the molecul o 34 i 35 , the transmission coeffi- 
cients can change dramatically, which modifies the slope 
near the Fermi level. Predicted thermopowers can also 
differ from experimental values even when the contact 
configuration is correct, due to approximations built into 
DFT. 

In contrast to BPDT, the CSW-479 molecule produces 
large magnitudes for the thermopower when the Fano 
resonance is within k^T of the Fermi energy. At low 
temperatures (T < 50if), the largest value of S is found 
to be approximately 200 /iVK" 1 and at higher temper- 
atures of approximately 670 /iVK" 1 . Furthermore the 
thermopower becomes negative as the angle increases and 
the Fano resonance approaches the Fermi energy from 
above. This suggests that thermoelectric properties can 
be tuned by modifying side groups and their coupling to 



the backbone. In this case it is also possible obtain a 
negative thermopower at high temperatures even though 
the Fermi energy is still inside the HOMO-LUMO gap. 
The sign change in the thermopower arises from a combi- 
nation of the asymmetric term Li and the intrinsic asym- 
metric shape of the Fano resonance. Therefore, the idea 
that the hole or electron-like character of conduction is 
determined by the sign of the thermopower is ambigu- 
ous for systems with Fano resonances close to the Fermi 
energy. 

As a demonstration of the generic nature of these pre- 
dictions, we repeated the above calculations for molecules 
connected to a gold atom on top of a gold electrode. The 
results were qualitatively similar, even though the details 
changed, such as the position of the Fermi energy, which 
shifted slightly closer to the HOMO level due to the re- 
duction of the charge transfer in the top configuration. 

Fig. (j4j) shows the calculated dimensionless figures of 
merit for a) CSW-479 and b) BPDT. Due to the low 
thermopower and conductance near the Fermi level the 
figure of merit for BPDT is low for all rotation angles, 
and therefore this molecule has a low efficiency of ther- 
moelectric heat transfer. This is not the case for CSW- 
479, where even at low temperatures a figure of merit 
ZT » 1 can be achieved for rotation angles where the 
Fano resonance is close to the Fermi energy. This demon- 
strates that the temperature difference can be mechani- 
cally tuned by adjusting the coupling between the local- 
ized states and the backbone. Such large figures of merit 
have been predicted for semi-conducting carbon nano- 
tubes 37 . However to date this is the largest predicted 
figure of merit value for an un-gated single molecular de- 
vice. 

We have shown there is a huge difference in the ther- 
moelectric properties of molecules with and without side 
groups. In contrast with electrical properties, the ther- 
moelectric properties of a molecular device depend on 
the shape of the transmission profile about the Fermi en- 
ergy rather than its absolute magnitude. By moving a 
Fano resonance through the Fermi energy, the shape of 
the transmission can be dramatically altered to produce 
both huge thermopower and figures of merit, correspond- 
ing to a high-efficiency thermoelectric cooling device. In 
practice, at high temperatures, this value could be lim- 
ited by the contribution to the thermal conductance from 
phonons 3 -^, which acts in parallel to the electronic con- 
tribution to the thermal conductance and is not included 
here. Electron transport through quasi-one-dimensional 
wires is also sensitive to disorde r 38 ! 39 , which could arise 
under ambient conditions, due to the presence of polar 
water molecules 19 . 
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